Have neutron stars a dark matter core? 
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Recent observational results for the masses and radii of some neutron stars are in contrast with 
typical observations and theoretical predictions for "normal" neutron stars. We propose that their 
unusual properties can be interpreted as the signature of a dark matter core inside them. This 
interpretation requires that the dark matter is made of some form of stable, long-living or in general 
non-annihilating particles, that can accumulate in the star [J|. In the proposed scenario all mass- 
radius measurements can be explained with one nuclear matter equation of state and a dark core 
of varying relative size. This hypothesis will be challenged by forthcoming observations and could 
eventually be a useful tool for the determination of dark matter. 
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Nowadays dark matter (DM) is an inevitable reality 
in both astrophysics [2[ and particle physics [3[ . The ev- 
idence for this non-baryonic ingredient of our Universe 
comes from observations at galactic and super-galactic 
scales, and it seems clear that it is the largely dominant 
fraction of matter. In this scenario it is to be expected 
that some DM should be present in stellar objects. De- 
pending on its nature, the mechanisms for accumulating 
DM are different, and it can result both from the stellar 
formation process and from subsequent accumulation via 
capture of DM particles during the whole stellar lifetime. 

On this basis, efforts have been made to investigate 
the possible observational consequences. To this aim, 
during the last years attention has been put to DM in 
relation to neutron stars (NS). The reasons are that NS 
should efficiently capture DM due to their high particle 
density, and new precise results from observations of NS 
are becoming more frequently available. Therefore, NS 
are considered potentially useful objects to help constrain 
DM models. In two recent papers jj,[5| the constraints on 
the properties of DM, in the form of weakly interacting 
massive particles (WIMPs), were revised by considering 
the effects on the luminosity of the accretion and self- 
annihilation of DM in NS. 

Here we follow another approach [l| and consider the 
potential effect on the masses and radii of NS equilib- 
rium configurations after accretion of DM. A key point is 
that the mass and radius of a NS are independent of the 
present DM accretion in the star. Instead, they depend 
on the whole DM capture process integrated over the 
stellar lifetime, and effects may be observable even if the 
star is not presently accreting DM. For some DM candi- 
dates there can be significant effects of DM on the equilib- 
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rium sequence of NS, i.e., on their masses and radii. The 
magnitude of this effect is proportional to the amount of 
DM trapped in the gravitational field of each NS, which 
depends on the individual history pj, starting from the 
formation of the progenitor star and continuing through 
the evolutionary phases until present age. Our hypothe- 
sis predicts significant effects on the mass and radius of 
some NS, in contrast to the small and presently unobserv- 
ablc modification of the luminosity expected from self- 
annihilation of WIMPs. Interestingly, there are now sev- 
eral observations of NS with mass-radius relations that 
could be in contrast with typical observations and predic- 
tions based on nuclear matter models |6H9|. Such obser- 
vations are used in attempts to constrain the equation of 
state (EOS) of NS matter [HJ. Here we consider observa- 
tional results for the four NS discussed in those papers. A 
combined analysis shows that there is a clear discrepancy 
between EXO 0748-676 [| and the group 4U 1608-52 j| , 
4U 1820-30 @ and EXO 1745-248 @. In fact, the author 
of Ref. [| claims that EXO 0748-676 rules out soft EOS 
for NS, while in the three more recent cases the same 
author with colleagues conclude that only a soft EOS is 
consistent with these mass-radius observations. In order 
to explain these data it has recently been suggested the 
existence of two populations of NS, compact and ultra- 
compact (composed partially or totally of quarks) [111 ] . 
see also fl2lll3J). 

An alternative explanation, which does not require 
novel exotic phases of dense matter, is that extraordi- 
nary compact NS have a dark matter core. This is not 
possible if DM is made of particles similar to the typical 
theoretical predictions for WIMPs, but only for stable 
or long-living particles that can accumulate in the centre 
of NS. A motivation and a detailed investigation of the 
effects of one such DM candidate on NS can be found 
in Ref. [lj and references therein. Here we repeat some 
of the key points. We consider stable DM particles that 
interact with ordinary particles by gravity only. Accord- 
ing to general relativity, all forms of energy are sources 



of gravity. In particular, the curvature of spacctime in- 
side a compact star depends on the energy-density and 
pressure distribution of matter. In presence of DM the 
metric is affected by both visible and dark components 
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In the Einstein field equations, the metric functions, \(r) 
and v(r), apply to both baryons and DM, and are mod- 
ified by the replacements ([T][2]). The standard equation 
for hydrostatic equilibrium separates for the two compo- 
nents 

d f = _[p( r ) + p(r)]fk (3) 

ar ar 

dr = -\pDM(r)+pDM{r)\-r, (4) 

because we assume that baryons and DM particles inter- 
act by gravity only, i.e., the gradient of the pressure of 
baryons does not exert a direct force on fluid elements 
of DM, and vice-versa. This result can be generalized to 
any number of fluids that interact by gravity [18[ . 

In a previous paper [l| we have solved these equations 
and studied the equilibrium properties of neutrons stars 
that have a core of mirror DM. This stable and self- 
interacting DM candidate emerges from the parity sym- 
metric extension of the Standard Model (SM) of particles 
[14| . In this model there exist a hidden set of particles 
that have the same physical properties as SM particles, 
but with "right-handed" weak interactions. These "mir- 
ror" particles interact with baryons essentially by gravity. 
For a review see Ref. [15| and for cosmological conse- 
quences Refs. [Ill [lTJ]- Since the microphysics of mirror 
DM is the same as that of baryons, we use the same 
nuclear matter EOS for the ordinary and dark matter 
components of the NS. In this special case, there are in 
fact two compact objects in the same location of space- 
time. One is visible and made of ordinary baryons, and 
the other is hidden and made of mirror baryons. In this 
scenario we find that the mass-radius relation is signifi- 
cantly modified in the presence of DM. This means that 
the equilibrium sequence of NS is non-unique and his- 
tory dependent, since it depends on the relative amount 
of baryons and DM. Interestingly, the mentioned recent 
observational results for the masses and radii of NS in- 
dicate that there might not exist a unique equilibrium 
sequence. 

In Figurc[T]wc plot the mass-radius relations for a mod- 
erately stiff EOS (for details see the original paper [l|) 
and a varying relative amount of mirror DM, tudm / ^tot, 
ranging from 0% to 70%. Values larger than 50% mean 
that there is more DM than ordinary matter, and that 
NS are the cores of bigger dark compact objects. Such 
objects should be extraordinary compact, making even- 
tual observations difficult to explain with other methods. 
A feature of this model is the shift of the equilibrium se- 
quence towards lower radii for increasing amounts of DM. 
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FIG. 1: Mass-radius relations for different dark matter (DM) 
contributions to the total neutron star (NS) mass, computed 
for the case of mirror DM. The numbers near the lines are 
mass ratios in percent, i.e., the mass of the DM core over the 
total stellar mass. Numbers larger than 50% are related to 
structures where the NS is the core of a bigger compact object 
made of DM. Observational results for four NS with 1-cr errors 
are included for comparison (see the text for references). 



Comparing with current observations, it appears that the 
less compact NS, EXO 0748-676, is well interpreted with 
models without DM, while for the other objects consid- 
ered here we need DM cores of at least 25% relative mass. 
In this way one can explain all the observational results 
with one nuclear matter EOS and varying amounts of 
trapped DM. If this picture is correct, it means that, 
under some conditions, DM can be efficiently captured 
by NS, or is present in significant amounts during the 
formation of some progenitor stars. In this scenario the 
opposite could also happen, i.e., that DM objects capture 
baryons. 

These results are calculated for mirror DM, but are 
qualitatively valid for other kinds of DM that could form 
stable cores inside NS. The basic prerequisites are that 
the DM particles are long-lived fermions with masses that 
are not too low or high compared to the baryon mass, and 
that they accumulate in sufficient number. Obviously the 
details of the effect on the equilibrium sequences depend 
on the particular DM model used, but the trend is general 
and leads to more compact objects in presence of DM. 
For example, it could be possible that, with some DM 
candidate, the observations are interpreted with a differ- 
ent EOS and smaller amounts of DM. Anyway, whatever 
the DM candidate that can accumulate inside NS, the ef- 
fect of adding a DM core is a softening of the EOS, thus 
mimicking that of other exotica as, e.g., quark matter. 

For example, DM particles with different masses and 
interactions than those considered here would affect the 
equilibrium sequence for a certain EOS in a different way. 
In general, whatever the DM candidate is, if it accumu- 
lates inside NS the effect is an apparent softening of the 
EOS, thus mimicking that of other exotica, e.g., quark 



matter. 

There is, however, a possibility to distinguish between 
this and other scenarios, since the amount of DM in NS is 
expected to vary and depend on the whole history of the 
star, especially on the environments from which it orig- 
inates and in which it lives. This will result in a spread 
of mass-radius measurements that cannot be interpreted 
with a unique equilibrium sequence. It seems unlikely 
that novel phases of ordinary matter can explain such an 
effect. 

Should our hypothesis be confirmed by future observa- 
tions of NS, a new possibility in DM research will arise: it 
will be possible to constrain DM properties by the inho- 
mogeneities derived using the estimates of DM inside NS. 
In fact, the comparison of masses and radii, derived from 
the observations, with the models computed for different 
DM contents, let us estimate the amount of DM in the 
core of each NS. This is the integrated result of the DM 
present at the formation of the progenitor star and the 
DM accreted during the whole life of the star (progeni- 
tor and NS). Correlating these results with the positions 
and orbits of NS, it should be possible to obtain informa- 



tion on the distribution of DM in the galaxy. Further- 
more, a reliable DM candidate should then also explain 
how it can accumulate in sufficient number in NS. Con- 
sidering the present estimated density for galactic DM, 
it's clear that a large amount of it cannot be accumu- 
lated only during the lifetime of the star, if we consider 
a smooth spatial distribution of DM. One or both of the 
following requirements should be verified: either the DM 
is present already during the formation process; or the 
density distribution of DM is highly non- homogeneous, 
with the possibility of events like mergers with compact 
astrophysical objects with stellar sizes made of DM. 

In conclusion, the recent observational results, which 
indicate that some neutron stars are extraordinary com- 
pact, can be explained with a nuclear matter equation 
of state if one allows for the possibility that dark mat- 
ter contributes significantly to the mass of some neutron 
stars. Improved and additional observational results will 
provide a crucial test of our hypothesis. If confirmed, 
observational results for the masses and radii of neutron 
stars will become a powerful tool for the determination 
of dark matter. 
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